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Abstract 

Thermal stratification and self- 
pressurization of partially filled liquid hydro- 
gen (LH 2 ) storage tanks under microgravity 
condition is theoretically studied. The tank is 
spherical and is subjected to a uniform and 
constant wall heat flux. It is assumed that a 
vapor bubble is located in the tank center 
such that the liquid-vapor interface and tank 
wall form two concentric spheres. This vapor 
bubble represents an idealized configuration of 
a wetting fluid in microgravity conditions. 
Dimensionless mass and energy conservation 
equations for both vapor and liquid regions 
are numerically solved. Coordinate transfor- 
mation is used to capture the interface loca- 
tion which changes due to liquid thermal 
expansion, vapor compression, and mass 
transfer at liquid-vapor interface. The effects 
of tank size, liquid fill level, and wall heat flux 
on the pressure rise and thermal stratification 
are investigated. Liquid thermal expansion 
tends to cause vapor condensation and wall 
heat flux tends to cause liquid evaporation at 
the interface. The combined effects determine 
the direction of mass transfer at the interface. 
Liquid superheat increases with increasing 
wall heat flux and liquid fill level and 
approaches an asymptotic value. A simplified 
correlation is developed for the maximum wall 
superheat as a function of system parameters 
using the numerical results. 


Nomenclature 


a 

acceleration of the tank 

Bo 

Bond number, par^/<7 

C P 

specific heat of liquid at constant 
pressure 

* 

C P 

dimensionless specific heat of liq- 
uid, C /C_ n 

^v.v’^p.v 

specific heats of vapor 

h f g 

latent heat of vaporization 

Ja 

dimensionless parameter, 

WS 

k 

thermal conductivity of liquid 

Lh 

dimensionless parameter, 
Po/(*o C p,o T o) 

Ls 

dimensionless parameter, 
2<r o/( r oPo) 

m e 

net mass transfer rate at the 
interface 

Nh 

dimensionless wall heat flux 
parameter, q w r 0 /(k o T o ) 

P,P 

tank (vapor) pressure, liquid 
pressure 

P* 

dimensionless tank pressure, p/p 0 

q w 

wall heat flux 


* 
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Qy vapor heating rate through vapor- 

wall boundary 

R* vapor bubble to tank radius ratio, 

r i/ r 0 

R Q value of R* at t* = 0 

r radius 

r i> r o vapor bubble radius, tank radius 

r* dimensionless radius, r/r Q 

T liquid temperature 

T* dimensionless temperature, T/T 0 

t time 

t* dimensionless time, t/(r 2 /a Q ) 

V y volume of vapor region 

Greek Symbols 

a liquid thermal diffusivity, k/(pC p ) 

a', fi f coefficients defined in Eq. (15) 

coefficients of liquid and vapor 
thermal expansion 

e*,r* transformed coordinates defined 

by Eqs. (14) and (15) 

K y isothermal compressibility of 

vapor 

p,p y liquid density, vapor density 

p* dimensionless liquid density, p/p Q 

a surface tension 

Subscripts 

o liquid property evaluated at the 

initial condition 

s evaluated at liquid-vapor interface 

v evaluated in vapor region 

w evaluated at tank wall 

Introduction 

The thermal environment in space can 
result in heat transfer through the thermal 
control system and increase the pressure of a 


cryogenic storage system. It is required that 
the tank be always kept within a speci- 
fied operating pressure during the entire 
space mission. Thus, knowledge of self- 
pressurization of a cryogenic fluid storage 
system is of great importance in the planning 
of space missions. Generally, the homoge- 
neous model, which assumes a uniform tem- 
perature in liquid and vapor phases, is used as 
a baseline to predict the tank pressure rise 
rate. An analysis of pressure control based on 
the homogeneous model has been performed 
by Lin, et al. 1 Ground experiments with 
liquid hydrogen (LH 2 ) have shown that the 
actual pressure rise rate spans a range from 
about l 2 to more than 10 times 3 faster than 
the homogeneous model predictions for a well 
insulated tank. This is due to thermal strati- 
fication in the liquid region which results in a 
subcooled bulk liquid and due to evaporation 
at interface. 

In a microgravity environment, the 
liquid-vapor interface is principally deter- 
mined by the surface tension force. The inter- 
face configuration also depends on the fluid 
properties, tank geometry, and liquid fill level. 
Usually, Bond number (Bo) and contact angle 
are used to characterize the interface configu- 
ration. Figure 1(a) shows the representative 
liquid-vapor configuration of a wetting fluid in 
a spherical tank for Bo < 1. The interface 
configuration on the ground where Bo > > 1 
is shown in Fig. 1(b) for comparison. For the 
tank with fully wetted wall, all the heat can 
be transferred from the tank wall to only the 
liquid region which is then superheated. Pres- 
sure increase is due to liquid thermal expan- 
sion and mass transfer at the liquid-vapor 
interface. This paper presents a theoretical 
solution of thermal stratification and self- 
pressurization of a spherical LH 2 storage tank 
in a microgravity environment. Liquid is 
treated as a compressible fluid such that the 
effect of thermal expansion on the solution is 
taken into account. The effects of tank size, 
liquid fill level, and wall heat flux on the 
pressure rise and liquid superheat are 
investigated. 
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Problem Formulation 


The physical system considered and coor- 
dinates are shown in Fig. 1(a). A spherical 
tank of radius r Q is partially-filled with a 
liquid. A vapor bubble of radius t- x is as- 
sumed to be located at the center of the tank 
such that the tank wall and the liquid-vapor 
interface form two concentric spheres. The 
vapor bubble configuration is representative to 
that of a wetting fluid in microgravity condi- 
tion for Bo < 1. The fluid in both phases is 
initially at rest and at a uniform saturated 
temperature T Q corresponding to the vapor 
pressure p Q . The tank wall is subjected to a 
uniform and constant heat flux q w . 

Vapor Region 


The mass transfer is related to the liquid 
temperature gradient at the interface as 
follows: 



( 3 ) 


Thus, Eq. (2) becomes 



ar 

dr . 


dr i 

~ Py 

dt 


r i^v*v(^v,v/^p,v) 


( 4 ) 


The vapor region is assumed to be at a 
uniform temperature. The pressure change 
rate of this homogeneous vapor region can be 
expressed as 4 


Clearly, the calculation in the vapor region is 
coupled with the liquid region through the 
liquid temperature gradient (5T/dr) 8 and 
through the change of interface location 
drj/dt. 


dp _ (/MVC PjV )+m e -p v (dV v /dt) 


dt 


vvV v « v (c ViV /c P)V ) 


( 1 ) 


where m e is the net mass transfer rate at the 
liquid-vapor interface: 


m 


e 


evaporation rate if m € > 0 
condensation rate if m e < 0 


In the present problem, there is no direct con- 
tact between vapor region and tank wall so 
that Q y = 0. Thus, Eq. (1) becomes 


Liquid Region 

The liquid is treated as a compressible 
fluid and the following assumptions are made: 

(1) The thermal conductivity, k, and the 
latent heat of vaporization, h fg , are assumed 
constant and fixed at the values corresponding 
to its initial condition. The variation of k 
and h fg in the range of temperature and 
pressure of interest is negligible. 

(2) All other liquid properties such as />, 
C p , /?, and a are considered variables. 

(3) The convective terms in the energy 
equation are neglected. 


dp = m e - ?v( dV v/ dt ) 
^ " ?vV v * v (C V)V /C piV ) 


( 2 ) 


(4) Energy dissipation is neglected. 


Thus, the energy equation in the liquid region 
is expressed as 


/> C P 


6T 

dT 


k o d r 2 dT 
r 2 dr( dr , 


+ 0T — 

dt 


( 5 ) 


Equation (2) shows that the pressure change 
results from the mass transfer at the liquid- 
vapor interface and the vapor volume change. 


3 



=0 ( 11 ) 


It should be noted that the liquid pressure P 
is different from the pressure p in the vapor 
region because of the surface tension effect. 
The relation between liquid and vapor pres- 
sure is given by 


f 1 . r*V dr* + Ja R * 2 


✓ R 

[dr*[ 


d 

dt* 


where the nondimensional variables are 
defined in the nomenclature. 


p= P -^ 

T: 


By using Eq. (6), Eq. (5) becomes 


„ 3T _ k o d 


20T 

r aTj 


( 6 ) 


+ /3T — 

dt 


2a 

P - — 

T: 


(7) 


Initial and Boundary Conditions 

The tank fluid is assumed to be initially 
at rest at a saturated pressure p 0 corre- 
sponding to an interface temperature T Q . In 
the present study, the liquid temperature is 
assumed to be uniform and equal to T Q . 
Thus, the initial conditions are: 


The rate of change of liquid mass is related to 
the mass transfer rate at the interface through 
the following conservation equation: 


dt 


f 1 ® 4jrr 2 a dr + 4jrr? 


( 

0T 



dr 

V / 


= 0 ( 8 ) 


T * = 1 
(p* = l 


att*=0 and R* < r* < I 


(12) 


The tank wall is subjected to a uniform and 
constant heat flux q w . The boundary condi- 
tion at the tank wall is 


Dimensionless F ormulation 

Dimensionless forms of Eqs. (4), (7), 
and (8) can be expressed as follows: 

Vapor : 


do* (3T fdr\ 

- 3 Ja(C ,/C J,,)\ _ 

df R*(p„«,) 


- 3(C p ,/C ) (dR’/dt’) 1 

R* (Po^v) 


(9) 


ar* 


\ dr* 


<lw r 0 

k T 
K o ± o 


= Nh 


(13) 


at r* = 1 for t* > 0 


It is noted that the temperature, T s (t*), is 
equal to the saturation temperature corre- 
sponding to the instantaneous vapor pressure. 
The instantaneous location of the interface, 
R*(t*), and its temperature, T g (t*), must 
satisfy Eqs. (9) and (11) simultaneously. 


Coordinate Transformation 


Liquid : 

,*c*?L* = J_ ± 

P * +2 -v.* 

r or 


dt * 


.2aT + 


dr* 


Lh(y9T) _ 

dt* 


♦ Ls 
P - — 


R*) 


( 10 ) 


Near the tank wall and the liquid-vapor 
interface, temperature gradients are expected 
to be large. Grid points should be clustered in 
these regions when the problem is solved 
numerically. Transformations are used which 
lead to a uniformly spaced grid in the compu- 
tational plane while points in physical plane 
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are unequally spaced. Also, the interface 
location is moving due to the liquid thermal 
expansion, vapor compression, and mass 
transfer at the interface. The fixed number of 
grid points on the interior are scaled to the 
motion of the boundary to provide adequate 
resolution. In order to accomplish the 
above two purposes, the following formula- 
tions are proposed by using the Roberts 
transformation 5 : 


Liquid : 

,-c 'El 

r P 


dr* 

♦ r ,*ae*ar* 


> 

1 de * 

d 

' *2 de *' 

dT * 

_ *2 t.* 

dr J 

de* 


de *. 


-p'c 


at* de* 


+ Lh(/?T) _ 
dr* 


* Ls 
P " — 

R* 

o 


( 17 ) 


r* = t* 


(14) 



-j 

1 T n, *2 

r p de +JaR 

r > 

dT * 

de*' 

e* = e*(r*,t*) = a' + (1 - a') 
r 

dr* J 

1 

0 

[de*i 

,U*J 


= 0 


In 


x 


ft' + (2a' + 1) (r* - R*)/(l - R*) - 2a' 


ft' - {2a' + 1) (r * - R*)/(l - R*) + 2a' \ 


In 


(£' + i) 


(A'-i) 


(15) 


where R* is a function of t* and /?' > 1. 
Equations (14) and (15) will transform an 
unequally-spaced physical plane (r*,t*) into a 
uniformly spaced grid in the computational 
plane (e*,r*). If a f = 0 the mesh will be 
refined near r* = 1 only, whereas, if 
os' = 1/2 the mesh will be refined equally 
near r* = 1 and r* = R*. If the stretching 
parameter is closer to 1, the meshes are 
more clustered. In the present study, the 
values of a ' and /?' used are 0.5 and 1.02, 
respectively. 

Based on this transformation, Eqs. (9) 
to (11) become 


Vapor: 


dP* _ O 


ar 




= 3 Ja 

dr* c v,v Py R *(p 0 /c v ) 


'de*' 


lar* 


/8 


(16) 


- 3 


P.V 




dR' 

dt* 


1 


R* Po K v 


(18) 


where 

r* = (1 ~ R*) 


{P+2a') 


fi' + l 




U'-iJ 


-ft' +2a' 







(2a' + 1) 

1 + 

' \ 

£' + 1 

(e *-«')/(! -a') 



U'-iJ 




+ R* 


(19) 


2ft'(l - a')[2a' + 1) 


— 

dt* 


In 

r 

ft' + 1 



[y - iJ 



f>' 2 ~ 


V - R*' 


[ 1 - R* 


(2a' + 1) - 2a' 



( \ 

- 

dR 1(1 - r*) 


dt* J 

(i - R*) 2 


(20) 


5 



2£'(1 - a ')( 2 a ' + 1 ) 


de* 

dr* 


In 


/?' + 1 
J' - l 


x 


A ' 2 - 
1 

1 - R* 


r * - R* 
1 - R* 


(2a' + 1) - 2a' 


( 21 ) 

Also, the boundary condition, Eq. (13), trans- 
forms to 


dT 


Nh 




U' ) 


dr* , 


Solution Algorithm 


( 22 ) 


(c) Check Eq. (18) for liquid mass con- 
servation. If satisfied, go to 
Step (2). Otherwise, adjust the 
guessed value of T g (r* n+1 ) and 
repeat Steps 1(b) and (c). 

(2) Solve Eq. (16) for vapor pressure 

p *( r * n +l) 

(3) Check if the difference between the vapor 
pressure p*(r* n+1 ) obtained in Step (2) 
and the saturation pressure correspond- 
ing to the interface temperature 
T g (r* n+1 ) obtained in Steps 1(a) to (c) is 
less than 10’ 4 . 

(a) If yes, the current solutions are the 
solutions for r* = r * n +l Set 
n = n+1 and go to Step (1) to 
start the calculations for the next 
time step. 


The governing equations (Eqs. (16) 
to (18)) are placed in finite-difference forms, 
along with the initial and boundary condi- 
tions. The Crank- Nicolson implicit scheme is 
chosen as the finite-difference method of solu- 
tion. The thermodynamic and transport prop- 
erties of the fluid are obtained by using the 
GASP code. 6 

An iterative procedure is necessary to 
determine the required location of the inter- 
face R*(t*) and the interface temperature 
T*(t*). The required R* and T g should 
satisfy that the calculated tank pressure based 
on Eqs. (16) to (18) is equal to the saturation 
pressure corresponding to T g . Assuming the 
solutions of t * = r* n are known (where n 
represents the n fch time step), the general 
procedure of carrying out the numerical com- 
putations for r* = r* n_hl is listed below: 

(1) For a given guessed value of R*(r* n+1 ), 

(a) Guess a value for the interface tem- 
perature T g ( r* n+ 1 ) . 

(b) Solve Eq. (17) for liquid tempera- 
ture distribution. 


(b) If no, adjust the guessed value of 
R*(r* n+ 1) and repeat the Steps (1) 
to (3). 

Results and Discussion 

The mathematical modeling described in 
the previous section shows that the solution of 
the problem depends on the fluid type and the 
parameters: t*, Ja, Lh, Ls, R Q , and Nh. 

For a given fluid and a given initial condi- 
tions, the parameters Ja and Lh are speci- 
fied. The pressure and temperature changes 
in the tank at any time t* are controlled 
only by the parameters Ls, R 0 , and Nh, 
which characterize the effects of tank size (r Q ), 
liquid fill level (ri/ r 0 ), and wall heat flux (q w ). 
For liquid hydrogen system of practical inter- 
est, the parameter Ls, is less than Iff 6 . 
Therefore, the term containing Ls in com- 
parison to other terms in Eq. (17) is 
neglected. In the present study, a liquid 
hydrogen system initially in a homogeneous 
state of 1 atm is considered. The ranges of 
R Q and Nh are selected to be 

0.4 < R* < 0.9 and 0.05 < Nh < 0.5 
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Figure 2 shows that, for a given R* the 
dimensionless pressure rise is higher for 
increasing values of Nh. Thus, for a given 
tank size and liquid fill level, tank pressure 
increases faster for higher wall heat flux q w . 
The pressure rise is due to the combined 
effects of the thermal expansion of the liquid 
and the mass transfer at the interface. The 
pressure rise rate as predicted by using a 
homogeneous model is also included in Fig. 2 
for comparison. The results show that tank 
pressure rise rate can be slower than the ho- 
mogeneous rate if there is no direct contact 
between tank wall and vapor region and heat 
is transferred to the liquid region only. This 
is due to significant liquid superheating. 

Figure 3 shows the effect of liquid fill 
level on the tank pressure rise. For a given 
Nh, tank pressure increases faster for lower fill 
level (higher R Q ) than higher fill level (lower 
R 0 ) except for the initial time period. Before 
the effect of wall heat flux reaches the inter- 
face, tank pressure change is basically gov- 
erned by the volume change of vapor region. 
Therefore, the equation of state yields 


Ap 


oc 


P 



(23) 


Since vapor volume change (-AV y ) due to liq- 
uid thermal expansion is approximately the 
same for both low and high fill levels, Eq. (23) 
shows that for a given tank the effect of liquid 
thermal expansion on the tank pressure rise 
should be more significant for higher liquid fill 
level. This liquid thermal expansion causes a 
higher pressure rise for higher fill level during 
the initial time period. Thermal expansion 
not only causes vapor compression to increase 
vapor pressure but also causes negative tem- 
perature gradient at the interface to induce 
interface condensation. Thus, liquid tempera- 
ture near the interface region increases. This 
impedes heat transfer from the wall to the 
interface to cause evaporation. Therefore, 
after the initial time period in which thermal 
expansion effect is dominant, the pressure rise 
becomes faster for lower fill level (as shown in 


Fig. 3) because heat has been transferred to 
the interface to cause evaporation. This phe- 
nomena can be better explained by using 
Fig. 4, which shows the effect of R q on the 
temperature gradient at the interface for a 
given Nh. For low fill level such as 
R Q = 0.9, the change in the interface tempera- 
ture due to the thermal expansion of liquid 
and compression of vapor is negligible. The 
temperature gradient at the interface is 
always positive and results in evaporation. 

On the other hand, for high fill level such as 
R 0 = 0.4, the pressure increases mainly due to 
the liquid thermal expansion and vapor com- 
pression. The interface temperature increases 
because of the increase in pressure, resulting 
in a negative temperature gradient near the 
interface and vapor condensation. It is noted 
that the oscillations of the data in Fig. 4 is 
probably related to the accuracy in calcula- 
tions of fluid and thermodynamic properties in 
GASP, which basically uses a curve-fitting 
technique. 

For low liquid fill level and low wall heat 
flux, the liquid compressibility and fluid prop- 
erty variations may be neglected. Under these 
conditions, the temperature distribution in liq- 
uid region can be predicted approximately by 
using the theoretical solution for transient 
conduction in a semi-infinite solid with con- 
stant surface heat flux, as long as the interface 
temperature has not responded to the effect of 
wall heat flux. 

4c 

For a given R q and Nh, liquid super- 
heat increases with time. Figures 5 and 6 
show that^the wall superheat AT*, where 
AT* = T w - T g , approaches an asymptotic 
value for a given R q and Nh. This asymp- 
totic value of wall superheat, (AT*) max , 
increases with increasing Nh and decreasing 
R q . A simplified correlation which is devel- 
oped for maximum wall superheat as a func- 
tion of R q and Nh using numerical results 
is given as follows: 


(at *)„ 


= 0.55 1 - R 


, 0.89 


Nh 


(24) 
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for (1 - RjNh < 0.15. For values of 
(1 - R Q )Nh greater than 0.15, Eq. (24) over- 
predicts (AT*) max . It should be noted that, 
although pressure rise rate is slower for a tank 
with higher liquid fill level (smaller R q ), the 
maximum liquid superheat (AT*) max is 
greater. A superheated liquid represents a 
metastable condition and may lead to a pres- 
sure spike either due to boiling or flashing of 
liquid. Therefore, Eq. (24) may provide useful 
information for the system designer to avoid 
the occurrence of pressure spikes due to liquid 
boiling, once the incipient boiling temperature 
is known. 

Conclusion 

An analysis of thermal stratification and 
self-pressurization of a spherical liquid hydro- 
gen tank with a spherical vapor bubble 
located at the center has been performed. The 
liquid-vapor interface configuration considered 
corresponds to that of a wetting liquid in a 
microgravity condition, i.e., Bo < 1. Based 
on the assumptions and conditions considered, 
the present study draws the following 
conclusions: 

- For a tank with fully-wetted wall, the 
pressure rise rate can be lower than the homo- 
geneous pressure rise rate because of superheat 
of the liquid region. 

- For a given liquid fill level, dimension- 
less tank pressure rise is higher for increasing 
values of Nh. Thus, for a given tank size and 
fill level, higher wall heat flux results in a 
higher pressure increase. 

- Tank pressure increases faster initially 
for a higher fill level because of the significant 
effect of liquid thermal expansion. However, 
after a period of time, the pressure increase 
becomes faster for a lower fill level tank 
because of the interface evaporation. 

- Liquid thermal expansion tends to 
cause vapor condensation and wall heat flux 
tends to cause liquid evaporation at the inter- 
face. The combined effects determine the 
direction of mass transfer at the interface. 


- Liquid superheat increases with 
increasing wall heat flux and liquid fill level 
and approaches an asymptotic value. For a 
tank with (1 - R 0 )Nh < 0.15, a simplified cor- 
relation, i.e., Eq. (24), is developed for the 
maximum wall superheat using the numerical 
results. 
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Figure 3. — Pressure rise as a function of time for 
various fill levels. 
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Figure 5. — Wall superheat as a function of time for various 
heat flux levels. 
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Figure 2. — Pressure rise as a function of time for various heat 
flux levels. 
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Figure 4. — Liquid temperature gradient at interface as a function 
of time for various fill levels. 
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Figure 6. — Wall superheat as a function of time for various 
fill levels. 


9 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 


Self-Pressurization of a Spherical Liquid Hydrogen Storage 
Tank in a Microgravity Environment 


AUTHOR(S) 

C.S. Lin and M.M. Hasan 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


WU- 593-21 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E -6759 


9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 105372 
AIAA -92-0363 


11. SUPPLEMENTARY NOTES 

Prepared for the 30th Aerospace Sciences Meeting and Exhibit sponsored by the American Institute of Aeronautics 
and Astronautics, Reno, Nevada, January 6-9, 1992. C.S. Lin, Analex Corporation, 3001 Aerospace Parkway, Brook 
Park, Ohio 44142; M.M. Hasan, NASA Lewis Research Center. Responsible person, C.S. Lin, (216) 977-0138. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified - Unlimited 
Subject Category 34 


13. ABSTRACT (Maximum 200 words) 

Thermal stratification and self-pressurization of partially filled liquid hydrogen (LH 2 ) storage tanks under micro- 
gravity condition is theoretically studied. The tank is spherical and is subjected to a uniform and constant wall heat 
flux. It is assumed that a vapor bubble is located in the tank center such that the liquid-vapor interface and tank wall 
form two concentric spheres. This vapor bubble represents an idealized configuration of a wetting fluid in micro- 
gravity conditions. Dimensionless mass and energy conservation equations for both vapor and liquid regions are 
numerically solved. Coordinate transformation is used to capture the interface location which changes due to liquid 
thermal expansion, vapor compression, and mass transfer at liquid-vapor interface. The effects of tank size, liquid 
fill level, and wall heat flux on the pressure rise and thermal stratification are investigated. Liquid thermal expansion 
tends to cause vapor condensation and wall heat flux tends to cause liquid evaporation at the interface. The com- 
bined effects determine the direction of mass transfer at the interface. Liquid superheat increases with increasing 
wall heat flux and liquid fill level and approaches an asymptotic value. A simplified correlation is developed for the 
maximum wall superheat as a function of system parameters using the numerical results. 



17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

10 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 






















